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Abstract 

We investigate the dependence of the chiral condensate (qq) on the temperature 
and quark density using the soft- wall holographic model of QCD, adopting geometries 
with black holes at finite temperature and quark chemical potential \i. We find that, 
for fi below a critical value, increasing the temperature the condensate decreases and 
vanishes at a temperature T ~ 210 MeV (at /i = 0). An analogous behaviour is 
observed increasing the chemical potential at fixed temperature. These results agree 
with the findings obtained by other methods. We also comment on the robustness of 
the results if geometries not involving black holes are adopted at low temperature, and 
an Hawking-Page transition is implemented. 



The strong coupling regime of quantum Yang-Mills theories, where perturbative meth- 
ods are not helpful, represents a problem requiring continuous investigations. This is 
particularly true for QCD, so rich of important soft effects the full understanding of which 
is still lacking. The knowledge of the dynamics of the QCD vacuum would shed light on 
aspects of primary importance, such as color confinement and spontaneous chiral symmetry 
breaking. [1] 

For Nf massless quarks quantum chromodynamics shows a global U(Nf)L®U(Nf) R = 
SU(Nf) L ®SU(Nf) R ®U(l) B ®U(l)A symmetry, whose SU(N f ) L ®SU(N f ) R part, at zero 
temperature and vanishing quark chemical potential, is broken down to vector SU(Nf)y 
(the baryon U(1)b symmetry remains exact, the axial U(1)a is anomalous). The sponta- 
neous breaking of this chiral symmetry is at the origin of effects such as the presence of 
Goldstone bosons (for Nf=2 they are the pions) and the large mass splitting between, e.g., 
light vector and axial-vector mesons, and is signaled by a non-vanishing vacuum expecta- 
tion value of the qq = q R qL + qLQR bilinear quark operator. 
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If the temperature T or the quark density (or the quark chemical potential /i) are in- 
creased, QCD is supposed to undergo a transition to a phase in which the chiral SU(N 
SU(Nf)n symmetry is restored. This produces changes of the hadron properties, and has 
the consequence that the quark-antiquark condensate (qq) vanishes [2]. 

The behaviour of (qq) versus temperature and quark chemical potential is therefore 
of fundamental interest, and indeed lattice QCD and effective field theory analyses have 
been employed to determine such a dependence. In particular, results have been obtained 
by lattice QCD as far as the T dependence is concerned [3J, while in the case of the 
quark density dependence, the problem connected to the non-positive definite fermionic 
determinant in the numerical evaluation of the discretized QCD partition function has 
not been fully solved [4] . Effective theories of QCD have also been used to establish such 
dependences, with a particular role played by chiral perturbation theory |5J, and by the 
Nambu-Jona Lasinio (NJL) model together with its Polyakov-NJL extension [6]. 

A method to access the nonperturbative regime of QCD has been developed in the 
recent past years, and is inspired by the Anti de Sitter/Conformal Field Theory (AdS/CFT) 
or gauge/gravity correspondence conjecture [THS]. It is denoted as holographic QCD (or 
AdS/QCD), and assumes the existence of a duality between large- A^ c QCD and a higher 
dimensional weakly coupled semiclassical field theory formulated on a curved AdS-like 
spacetime. Using perturbative methods in the weakly coupled theory, one can obtain 
information about the dual strongly coupled QCD, since the correspondence implies that 
the partition functions of the two theories coincide, under suitable boundary conditions [9], 
and that the Minkowski spacetime on which the gauge theory is defined can be identified 
with the boundary of the AdS space. 

Our purpose is to study the dependence of the chiral condensate on the temperature and 
quark density in a holographic framework denoted as soft-wall model |10|, 11 1] (investigations 
in different holographic approaches can be found, e.g., in [121 [13]). The soft-wall model is 
defined in a non-dynamical AdSs spacetime; finite temperature and quark density effects 
are described considering geometries with black holes. In the case of finite temperature, 
the line element in the gravity model is chosen as the AdS-Schwarzschild (AdS/BH), 



R is the AdSs radius and z the fifth (holographic) coordinate (e < z, with e — > + ). Zh 
is the position along z of the horizon of a black hole, and it is related to the Hawking 
temperature T by the condition 




(1) 



with 




(2) 



T is identified with the temperature of the QCD boundary theory. Notice that we use 
capital latin indices for 5D coordinates, and Greek indices for 4D coordinates. With this 
choice of the metric we are assuming that the AdS/BH geometry is the stablest one, i.e. 
the lowest-action gravity configuration for all values of the temperature; we shall come to 
this point in the following. 

The soft- wall model involves a complex scalar (p = 0-form) field X(x,z) [TU], whose 
mass is fixed by the AdS/CFT prescription: m 2 R 2 = (A — p)(A + p — 4) [9]. This field is 
dual to the A = 3 dimension QCD qq operator. The gravity theory also includes 5D gauge 
fields, Lm and Rm, dual to the 4D A = 3 left-and right-handed currents j£ = q^l^lL and 

= QRl^QR associated to the global SU{Nf)i and SU{Nf)n symmetries, respectively: 
hence, the global chiral symmetry of the QCD boundary theory is gauged in the bulk. The 
5D meson-action defining the holographic model reads (we set to 1 the AdSs radius R) 
PI [H] 
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where g is the determinant of the metric ([!]) and Fl^r are the field strenghts F^ IN = 
pMNarpa = qM l n _ qN L M _ i [l m ,L n ], and an analogous expression for F R . The 

constants ky m and g\ have been fixed in [TU [15] : kyu = and g\ = - . The dilaton 



N c *° 4 

<j>(z) is a background field introduced to account for the breaking of conformal invariance 
in the IR, with the functional form 4>(z) = c 2 z 2 chosen to produce confinement [16J and 
linear Regge trajectories for light hadrons [10]. The scale parameter c is fixed from the 
mass of the p meson to the value c = m p /2 = 389 MeV [10] (in the following we set c = 1 
and express the dimensionful quantities in units of c). 

The field X is in the bifundamental representation of SU(Nf)i ® SU(Nf)n, while the 
gauge fields are in the adjoint representation of the corresponding Lie algebras: Lm = 
L M T£ and Rm = Rm^ri with T£ R the generators of SU(Nf)j J> ^\ the trace is taken over 
gauge indices. The covariant derivative acting on X is defined as TjmX = 8m X — iLmX + 
iXRm, and it is the way in which the two fields Lm and Rm are coupled to each other. 
Therefore, the scalar sector of the theory, i.e. X, is responsible of chiral symmetry breaking 
in the model (for X = chiral symmetry is restored) [T7], [18]. The fluctuations around 
the vacuum expectation value of X describe light scalar mesons |10} IT4] . To complete the 
model, a Chern-Simons term should be added to describe the chiral anomaly [19j; the dual 
of another A = 3 spin one operator = qa^q has also been proposed to be included in 
the gravity action [20]. Prom now on, we consider the case of two quark flavours, putting 
Nf = 2. The result can be easily generalised to other values of Nf. 

In the action Q the integration over the holographic coordinate z is extended up to 
the black- hole horizon z^, which represents an IR bound and sets the size of the bulk. We 



1 Lm and Rm are in the singlet representation of SU(Nf)n and SU(Nf)L, respectively. 
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assume that the behaviour of (qq) is due to the interaction with the geometry and the 
dilaton. This means that we only consider linearized equations for the field X, ignoring its 
coupling to the gauge fields, as well as other contributions coming from a possible potential 
term V(X) which could be added to the action Eq.Q [TO] . 

Defining the scalar field as X(x, z) = Xo(z)ljVf e* 7 ^' 2 -' (hence neglecting the fluctuations 
around the Xq configuration, while tt(x, z) = TT a (x, z)T a represents chiral fields), where ljy, 
is the Nf x Nf identity matrix, we can identify Xq(z) as the v.e.v. describing the dynamics 
of the condensate (qq)' as it can be inferred from dimensional analysis, the coefficient of 
the term proportional to z 3 in the expansion of Xq for z — > is indeed proportional to (qq) 
|21j.|15j. To determine such a coefficient, we work out the equation of motion for Xq(z) 
from the action Q: 

*S(z, %) " 2z2f{z) l\~ f{z) X' Q (z, z h ) + -^—,X (z, z h ) = (5) 

in which we have explicitely indicated the dependence of Xq on the position z^ of the 
horizon, hence on the temperature. The primes denote derivatives with respect to z. At 
T = (i.e. for Zh — > oo and f(z) = 1) this equation can be solved analytically and has 
only one regular solution in the (large z) IR region, 



X T =\z) = r ^pzU{^A^ (6) 

with U the Tricomi confluent hypergeometric function. The solution should also be mul- 
tiplied by a factor describing the scaling of Xq with N c |19| [22], which however is not 
relevant for the present analysis. The coefficient m q is identified with the quark mass by 
the condition that Xq/ z must tend to the source of the qq operator in QCD as z — > [9]. 
Eq. ([5]) also admits another linearly independent solution, which is singular for z — > 00; it 
corresponds to a divergent on-shell action Q and therefore it is discarded. 

In Eq.([6j the chiral condensate, the coefficient of the z 3 term in the expansion of X(z) at 
small z, and the quark mass, the coefficient of the linear z term, are proportional, a relation 
which does not hold in QCD. Indeed, expanding ([6]) for z->0we get the asymptotic form 

XT=\ Z ) ^m qZ + m q ( - l og (2) ) z 3 + m q z 3 log(z) + 0(z 5 ) . (7) 



- log(2)j z 3 + m q z 3 log(z) + 0(z 5 ) . 



This result is a drawback of the soft-wall model, and it does not allow to separately de- 
scribe explicit and spontaneous chiral symmetry breaking. In other holographic approaches, 
namely the hard- wall model [15] , the coefficients of the z and z 3 terms in the small z expan- 
sion of Xq are not related to each other; however, in those cases the chiral condensate is an 
external input and cannot be extracted from the dynamics of the gravity system (studies 
of temperature effects in the hard- wall holographic model can be found in [23, 24]). A way 
out in the soft-wall framework includes the possibility of adding a potential term to Q 
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|1Q[ 125] . together with other solutions [26]. In spite of this, it is possible to have hints on 
the dependence of (qq) on the temperature at T ^ 0. 

Near the black-hole horizon (z ~ Zh) Eq.Q admits two solutions, a logarithmically 
divergent solution and a regular solution which behaves as 



X {z,z h ) 1 - - 



z 

Zfl 
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6 I 



z 

Zh 



+ 0((z h -z) 3 ). (8) 



This regular solution for T — > continuously reproduces the zero temperature expression 
of Xq given by Eq. Q , therefore it can be used to fix the profile of Xq close to z^ . On the 
other hand, at low values of z, the two independent solutions have asymptotic form 



X i(z,z h ) 
X 02 (z,z h ) 



^> z + ^z 3 + z 3 log(z) + 0(z 5 ) 



> z 3 + 0(z 5 ), 



(9) 
(10) 



with the term proportional to z 3 log z having the same coefficient of the linear z term. The 
dependence on z^ (hence on the temperature T) does not enter the \ow-z behaviour up 
to 0(z 5 ). The general solution has then the form Xq = AXqi + B Xq2 and the condition 



z 

X is 



fixes the coefficient of (9) to A 



ni 



Therefore, the low-z dependence of 



X(z,z h ) > AX 01 (z,z h ) + B X 2[z,z h 

z 3 log(z) + (~ 
m q z + m q z 3 log (2) + az 3 



m q z + m q z 3 log(z) + ( - m q + B ] z 3 



(11) 



where a = | m g + -B is proportional to the chiral condensate (qq) through a temperature- 
independent coefficient; it is possible to find this relation from a matching condition with 
QCD of a one-point correlation function of a scalar operator [T3] (see appendix). Such a 
proportionality relation ensures a correct N c scaling of the field Xq (~ N®), the quark mass 
(~ N®) and the chiral condensate (~ N c ) [22]. In the soft-wall model, the coefficient cr(T) 
must be fixed through a suitable IR boundary condition. As such boundary condition, we 
require the general solution to match the behaviour pj as z — > z^, which ensures regularity 
and continuity with the T = case. To do this, we use a numerical shooting method, where 
the shooting parameter is the coefficient a in (11). The value of this parameter turns out 



to be different for different values of temperature, and applying the method for different 
T's we obtain the dependence of a (and, so, on the chiral condensate (qq)) on T in this 
model. m q is the bare quark mass and it is an external parameter in the soft-wall model 
[10] : we keep it fixed with respect to temperature, as e.g. in [1 2|. [13] and from this point 
on, we set m q = 1 for simplicity. 
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The resulting cr(T) is drawn in Fig. [T] Starting from the value at T = and increasing 
the temperature, after a slight increase to a local maximum at T/c ~ 0.31, the condensate 
drops to a = at T/c ~ 0.54. Numerically, this temperature corresponds to T ~ 210 MeV. 
This behaviour is similar to what is commonly expected: for example, in some finite tem- 

/ ( T' 

perature QCD sum rule analyses, the condensate profile o(T) = <t$ 



1 



is used 



(at least close to the critical temperature T c , with ao the value at T = 0) [27] . Analogous 
results are obtained by Schwinger-Dyson calculations [28] and in models of QCD [291 130] . 
Since for T ^ T c the chiral condensate vanishes, for massless quarks one has chiral sym- 
metry restoration and T c is the temperature at which the chiral phase transition occurs, 
with a a critical exponent. Lattice QCD analyses also find a decreasing chiral condensate, 
which vanishes at T ~ 180 — 200 MeV [3|. A summary of the results obtained through 
different lattice approaches can be found in |31j. 

In the soft-wall model chiral symmetry is never restored, since the quark mass remains 
different from zero; explicit symmetry breaking occurs and the chiral limit cannot be prop- 
erly investigated. On the other side, the fact that the behaviour of a for T ^ T is similar 
to the one expected in QCD for T ^ T c , is an interesting feature in view also of possi- 
ble modifications aimed at an improved description of chiral symmetry within this model. 
Extending T to higher values, a becomes negative, therefore T represents the maximum 
temperature below which the model can be suitably applied to hadrons. This temperature 
is higher than the one where melting is observed in the spectral functions of the lightest 
hadrons [51] . 

Let us now comment on the geometry assumed to be dual to finite temperature QCD. 
In frameworks in which a dynamical geometry is considered (in particular in the case of 
the gravity dual of N = 4 SYM), it has been observed that two different solutions of 
the Einstein equations in the bulk have to be taken to describe the finite temperature 



<r(T) 
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Figure 1: Dependence of a in Eq. ( 11 ) on the temperature T. Both quantities are expressed 
in units of the scale c. 
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field theory on the boundary, the AdS/BH and the thermal- AdS (th-AdS). The latter is 
characterized by AdS metric and compact periodic Euclidean time dimension, with size 
proportional to 1/T. Depending on the temperature, the two geometries correspond to 
different stationary points of the gravitational action, th-AdS being the stablest solution 
at low temperature, AdS/BH the stablest one at high temperature. A Hawking-Page 
transition between th-AdS and AdS/BH is interpreted as the gravity dual of the QCD 
deconfinement transition. 

In the soft-wall model, assuming that the th-AdS+dilaton and AdS/BH+dilaton ge- 
ometries correspond to exact solutions of gravity equations and considering the on-shell 
gravity action, the existence of a Hawking-Page transition has been envisaged [32], with 
a deconfinement temperature Tjpp ~ — g 47 ~ 0.49 If the th-AdS and the AdS/BH 
geometries were adopted in the present calculation of the chiral condensate, following the 
method discussed above, the chiral condensate would remain constant up to T = Thp, 
and then would follow the curve in Figjl] for higher temperatures. The small increase of 
the condensate with temperature shown in Figjl] (not observed on the lattice and in chiral 
perturbation theory) would be absent, and the profile for a at T = Thp would be dis- 
continuous. However, the assumption that the geometry and the dilaton field effectively 
result from the solution of a coupled gravity-dilaton equation is crucial for the existence of 
the Hawking-Page transition [32J. On the other hand, the soft-wall model adopted here is 
phenomenologically constructed, the dilaton profile being suitably chosen to reproduce em- 
pirically observed features of the dual theory, and its dynamical origin is not known. Even 
standing these limitations, the soft-wall model can be considered a minimal benchmark in 
the attempt to get information on important features of the QCD vacuum. 

To investigate the dependence of the chiral condensate on the quark density, we need to 
introduce the quark chemical potential in the holographic framework. Finite quark density 
effects can be described by adding to the action, in the QCD generating functional, the 
term f d A x /J,q*q, where q^q = qj°q is the quark number operator and [i the quark chemical 
potential. Therefore, [i can be considered as the boundary value of the time component 
of a £7(1) gauge field in the bulk theory. The geometry of the bulk coming from the 
interaction between a space-time with negative cosmological constant (the AdS space) and 
a £7(1) gauge field results in an AdS/Reissner-Nordstrom (AdS/RN) space, an AdS space 
containing a charged spherically symmetric black hole. For vanishing spatial components of 
the £7(1) gauge field Ai(z) = (i = 1, 2, 3, z), the small z expansion of the time component 
Aq(z) can be written as 



where q is the charge of the black-hole and k a dimensionless coefficient which can be 
considered a parameter of the model [35J. k scales as y/Nc, and its (model dependent [36J) 

2 A critical temperature Tc — 210 MeV, associated to the deconfinement transition, has been obtained 
in a soft-wall model with AdS/BH geometry and no Hawking-Page transition, analyzing the static quark- 
antiquark pair potential 33 . 




(12) 
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numerical value represents a parameter useful to phenomenologically implement subleading 
effects in the large N c limit. The results in the following refer to k = 1. The line element of 
the AdS /RN gravity theory has the same form as in , with the black-hole factor modified 
to 46 

/(*) = 1-(1 + Q 2 ) (±) +Q 2 [£) 0<Q 2 <2, (13) 

where Q is Q = qz\ [35J. The quantities Q and zy L are connected to the temperature T, 
defined by Q, and to the chemical potential \x: 

TTZh v 2 ; 

(J, = K— (14) 

with the second equation coming from the condition that ^4q vanishes at the horizon. 
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Figure 2: Chiral condensate versus the temperature T at different values of the quark 
chemical potential fi (a), and versus \x at different values of T (b). In (a) the curves 
correspond to fj,/c = 0.02 (plain, red), 0.8 (blue, dashed), 1.2 (green, dotted) and 1.6 
(black, dot-dashed). In (b) the curves correspond to T/c = 0.05 (plain, red), 0.3 (blue, 
dashed), 0.4 (green, dotted) and 0.5 (black, dot-dashed). The parameter k is set to k = 1. 

Following the same procedure outlined above to solve ([5]), we find that as z ~ 0, the 
two linearly independent solutions follow Eqs.(9fl~0 ). On the other hand, near the horizon, 
the regular solution now behaves as 



X(z,z h ) 



^ 1 



+- 



3 
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z h, 
->2 i Atr& 



2Q 2 + 4(Q 2 -2)z 2 h ) 
16(Q 2 - 4) 



z 

Zh 



+ 0{{z h -zf). (15) 
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We numerically solve the equation of motion (|5| tuning a in Eq. (11) to obtain the asymp- 
totic dependence (15) at z ~ Zh, and the results of the calculation are collected in F ig. [2 



The behaviour of a versus T for different values of chemical potential /x is plotted in 2(a) 



for all the considered values of the chemical potential, increasing the temperature the chira 
condensate drops to zero at values of the temperature that depend (decrease) on \x. An 



analogous behaviour is found exchanging T with /i, as shown in Fig. 2(b). These depen 



dences agree with the results obtained by different methods, for example in approaches 
based on the composite operator formalism [37]. The physical value of fi for which, at low 



temperatures, the chiral condensate vanishes depends on the parameter k in (12): using 
the value k = 1/2 fixed (in a different model) in [35J, it corresponds to fi ~ 350 MeV. 

Analogously to the fj, = case, also at finite density, if one considers a dynamical 
geometry and dilaton, it is necessary to establish which is the stablest geometry for each 
value of T and fx. This problem has been studied in top-down approaches [12] and in the 
bottom-up hard- wall model [36, 38J. In [12J it is shown that for fj, > the thermal-AdS 
solution is inconsistent, and the geometry always contains a black-hole. In particular, a first 
order phase transition between two black- hole geometries occurs at a critical T c (/j,), up to a 
certain value of \x. In the hard-wall model it is found that a transition between a "charged- 
th- AdS" solution and the AdS /RN one could take place [361 EH] • In our phenomenological 
model the geometry and the dilaton are fixed background configurations, the existence of 
transitions between different geometries being related to the unknown solutions of gravity 
equations: for this reason we only consider AdS/RN for all values of T and fi. 

In Fig. [3] we plot the values of the temperature f where a vanishes, versus the chemical 
potential \i. It is similar to the one expected in the case of transitions to a phase with 
restored chiral symmetry: in the case of the soft- wall, the curve depicted in Fig. [3] indicates 
the region of T and fi in which the model can be used to describe hadrons. 

To conclude, studying (qq) cx a at finite temperature and quark density using the 



0.7 r 
0.6 1 




Figure 3: Line in the T — \i plane where a vanishes. The parameter k is set to k = 1. 
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soft-wall holographic model of QCD and geometries with black holes, we have found a 
behaviour of the condensate versus the temperature in qualitative agreement with the 
results of other approaches. In particular, increasing the temperature cr(T), after a slight 
increase up to a maximum, drops to zero at a value T ~ 210 MeV which is close to other 
determinations. For finite quark chemical potential the behaviour is analogous, and the 
line in the temperature/chemical potential plane in which a vanishes is similar to the line 
envisaged for QCD. Such results encourage the hope that a dual theory of QCD can be 
formulated. 
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Appendix 

A proportionality between the coefficient cr(T) and the chiral condensate is found here, 
showing that all the temperature dependence of (qq) is contained in <r(T), while the co- 
efficient of proportionality is temperature-independent and ensures a correct scaling of 
quantities with the number of colours. The chiral condensate can be computed within 
holographic approaches through the relation 



(qq)= lim W— (4) . (16) 




being \Jg) the vacuum expectation value of the quark scalar current. In order to evaluate 
this quantity, let us consider the on-shell action for the scalar field 



i r „ p-'M 2 ) f( z ) 

- / (Tx ^^Tr[X(x,z)d z X [X, z ) 

k Y M J z 6 



(17) 



z=0 

where 

X(x,z) = {X {z)l + S(x,z))e 2i7T . (18) 

Xq(z) is the expectation value we have considered in this paper, ir is the chiral field, and 
S(x,z) is the fluctuation describing scalar mesons |14j : 

S(x, z) = J d A y S(x - y, z) S a (y) T a , (19) 

S(x — y,z) being the bulk-to-boundary propagator and Sq the source, according to the 
holographic dictionary [S]; T a (a = 1, ..8) are the Gell Mann matrices and T° = 1/ \J2nf. 
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Then in the holographic model the relation holds: 



5S 

5S$ 



(20) 



and, in the Fourier space, defining S(q, z) = J cfix e ig ' x S(x, z), the vev of the scalar current 

reads: 

y/2n~f f(z) 



2 kyu 



X (z) S'(0,z)+X' o (z)S(0,z) 



(21) 



2=0 



where the prime indicates a derivative with respect to z. At finite temperature and for low 
values of the fifth coordinate z, the scalar field S(q, z) can be written as: 



S(q,z) 



2-S>0 



> z + A(T) z 3 log z + B{T)z 3 + 0{z" 



(22) 



where we have used the boundary condition S(q,z)/z z ~ > °> 1, and A(T) and B(T) are 
coefficients we do not need to specify for the scope of the calculation. The low-z behaviour 
of Xq is 

(23) 



X (z) ^> m g z + m q z 3 log z + a{T)z 3 + 0(z b 



where only cr(T) depends on the temperature, as shown in the paper. Plugging (22) and 
(23) into (21), we find the proportionality relation between the coefficient cr(T) and the 
chiral condensate^] 

4 TV 
(qq) = - T —a(T)=--^a(T), (24) 



kyu 



provided cr(T) is finite in this limit. This result shows that the temperature dependence 
of the chiral condensate is only governed by c(T), the proportionality coefficient N c /(4tt 2 ) 
not depending on T. It is worth stressing that such a coefficient ensures the correct scaling 
of all the quantities involved in Eq. (23) with N c . This solves the problem of the scaling of 



the vev with N c pointed out in [22], since all the coefficients in the expansion (23) scale as 
N°, while (qq) scales as N c . 
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